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Introduction 


A project on droplet combustion, pursued jointly with F.L. Dryer of Princeton University, has 
now been in progress for many years. The project involves experiments on the burning of single 
droplets in various atmospheres, mainly at normal atmospheric pressure and below, performed in 
drop towers (e.g., ref. 1) and designed to be performed aboard space-based platforms such as the 
Space Shuttle or the Space Station. It also involves numerical computations on droplet burning, 
performed mainly at Princeton, and asymptotic analyses of droplet burning, performed mainly at 
UCSD. The focus of the studies rests primarily on time-dependent droplet-burning characteristics 
and on extinction phenomena. 

The presentation to be given here concerns the recent research on application of asymptotic 
methods to investigation of the flame structure and extinction of hydrocarbon droplets. These 
theoretical studies are investigating the extent to which combustion of higher hydrocarbons - 
heptane, in particular - can be described by four-step reduced chemistry of the kind that has 
achieved a good degree of success for methane flames. The studies (refs. 2, 3, and 4) have progressed 
to a point at which a number of definite conclusions can now be stated. These conclusions and the 
reasoning that led to them are outlined here. 

Asymptotic Descriptions of Diffusion Flames 

A mixture fraction Z (ref. 5), defined to be zero in the oxidizer and unity in the fuel, serves as 
an attractive independent variable for analyses of diffusion-flame structures because its adoption 
enables the problem to be formulated in a universal manner that does not depend on the specific 
geometrical configuration. However, a unique Z exists only if all relevant Lewis numbers are 
unity. It has been argued (ref. 3) that for droplet burning, in the transport zones away from the 
stoichiometric mixture fraction Z c , Lewis numbers of unity are reasonable approximations, mainly 
because those for oxygen and products in the outer zone are close enough to unity, while in the fuel 
zone only the fuel specific heat is relevant (since no other species are transported there), and this is 
large enough to counteract the small fuel diffusion coefficient, again giving the Lewis number close 
to unity. It is important to account for Lewis numbers different from unity in the thin reaction zone 
centered at Z = Z C1 but this can be done without difficulty and without affecting the definition 
of Z, so long as this zone is thin in the mixture-fraction coordinate. Thus, the droplet-burning 
problem is formulated in terms of the mixture fraction (ref. 2, 3, and 4). 

In the mixture-fraction formulation, the simplest asymptotic description is the Burke-Schumann 
reaction-sheet structure, shown by the straight lines labeled “equilibrium” in figure L This limit 
corresponds to an infinite Damkohler number for one-step, irreversible chemistry and provides 
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the outer, transport-zone solutions to which the reaction-zone solutions are matched. Figure 1 
illustrates that in the reaction zone there is a one-parameter family of solutions, the parameter 
being Xo the value at stoichiometry of the scalar dissipation rate 


X = 2D\VZ 


( 1 ) 


where D is the (local, variable, common) diffusion coefficient of the outer solutions. As x c increases, 
the peak temperature decreases, and eventually an abrupt extinction is approached, at a critical 
value of Xc- Figure 1 illustrates this qualitatively on the basis of activation-energy asymptotics 
(AEA); after extinction the solution is the frozen one given by the dashed lines. 

When the mixture-fraction formulation is employed, the distinction between different geomet- 
rical configurations arises only through the expression for the scalar dissipation, x- For a stagnant 
mixing layer of thickness, l 

X = 2 D/l, W 

while for constant-density counterflow mixing with a normal-component strain rate Ko m the 
oxidizer stream , 

x = (Ao/7r)/exp |2 [erfc- 1 (2Z c )j | ( 3 ) 

in which erfc" 1 denotes the inverse (not the reciprocal) of the complementary error function. It 
has recently been shown (ref. 6) that for small Z c , the usual case encountered in practice, variable- 
density effects can be included reasonably by placing on the right-hand side of equation (3) t e 

additional factor , , -i 

/ = |3 (\fpoJpc + l) 2 j / J4(2 \fpo/pc + 1)J * W 

where po is the density of the oxidizing stream and p c that at the stoichiometric flame sheet. For 
droplet burning (ref. 2) it is found that 
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ln( 1 + B) 


}■ 


( 5 ) 


where d t is the droplet diameter, and B is the transfer number for the droplet combustion process, 
given in an approximation of constant specific heat c p by 


B = {c^Too - T i) + QYo 2 oo/ [(3« + 1)WV2]} 


( 6 ) 


where T,W, and Y denote temperature, molecular weight, and mass fraction, the subscript oc 
identifies the ambient atmosphere, L is the energy required to vaporize a unit mass of fuel, Q 
stands for the heat released per mole of gaseous fuel in combustion, and n represents the number of 
carbon atoms in the alkane. By use of equation (5) at Z = Z c , the droplet diameter at extinction 
du, can be obtained from y„. the valu e of Xc at e xtinction, which has been calculated without 

reference to the specific geometrical configuration. 

Figure 2 iUustrates how figure 1 is modified by reduced chemistry. This figure pertains to a two- 
step approximation in which the reaction zone is split into two parts, a fuel-consumption zone and 
an oxidizer- consumption zone. The asymptotic analyses have addressed flame structures m two- 
step, three-step, and four-step approximations. Qualitative illustrations of resulting concentration 
and temperature profiles in physical space are shown in figure 3, in which C 2 H 2 is one representative 
of numerous intermediate fuel species. 
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Reduced Chemistry 


Systematic reductions to four-step mechanisms for alkanes result in the reduced chemistry 


C n H2n+2 + 7 H + \l3 - (n + 1 + ^7)] H 2 O + jn - 

(3 ~(n+ 1 + ^ 7 )j| 0 2 — ► olCO + @H 2 + (n - a)C0 2 , 


CO + H 2 0 s co 2 + h 2 , 
2 H + M -> H 2 + M, 
0 2 + 3ff 2 ~ 2H + 2 H 2 0. 


(7) 

( 8 ) 
(9) 

( 10 ) 


Equations (7), (8), (9), and (10) describe the fuel- consumption, water-gas shift (or CO-consumption), 
radical-recombination, and oxygen-consumption (radical- production) steps, respectively, and they 
proceed mainly at the rates of the elementary steps 


C n H 2n+2 + H - c n n 2n+l + h 2 , (11) 

CO + OH - C0 2 + II, (12) 

H + 0 2 + M -+ H0 2 + M, (13) 

and 

H + 0 2 -0H + 0, (14) 

respectively. The values of the constants a,/3, and 7 in equation (7) depend on the chemical-kinetic 
description of the fuel chemistry. Studies (refs. 3 and 4) have shown that theoretical predictions are 
especially sensitive to the value of 7, the number of radicals destroyed for each fuel molecule that 
is consumed. It may be expected that the larger 7 is, the lower will be the radical concentration 
and the less robust the flame. 

For methane (n = 1), the summation 


CH a + II 
CH 3 + 0 
CH 2 0 + H 
CHO + M 

H + OH 

H + II 2 0 
CIU + 2H + II 2 0 


- CH 3 + H 2 
-» CII 2 0 + II 

-> cho + ii 2 

-> CO + H + M 

- o + h 2 
-» oh + h 2 

-> CO + 4H 2 


with the steps dictated by the need to cancel all steady-state species through their fastest reactions 
results in a = 1,0 = 4, and 7 = 2. Even excluding the second of these steps by assuming 
(erroneously) that 2CH 3 + M — * C 2 He + M is very fast and that ethane chemistry therefore 
controls the fuel consumption gives only a = %,(3 = 3, and 7 = 1, that is, the number of radicals 
consumed is decreased by just a factor of two. Consequently there is relatively little significant 
uncertainty in step (7) for methane. However, for a larger alkane, such as heptane, there is a large 
difference, depending on whether the fuel chemistry follows a propene path after step (11), namely, 

CVffis — > CH 3 + 2C 3 IIq (15) 
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(16) 


(overall), as proposed by Warnatz, or an ethene path, 

CjHxs — ► CH 3 4* 3C 2 ^ 4 

(overall) as favored by Dryer. The propene path gives a = 7, (3 = 20, and 7 = 10, while the ethene 
path gives a = 4, (3 = 13, and 7 = 2, a factor of five difference in the number of radicals consumed. 
This large factor turns out to have a large influ enc e on the scalar dis sipation at extinction because 
the much greater extent of consumption of radicals by fuel in the prope ne pa th leads to extinction 
much earlier. Since uncertainties in the detailed fuel chemistry of higher hydrocarbons do not 
enable a selection to be made between the propene and ethene paths, for example, at present, it 
is concluded that further study of elementary steps and mechanisms in alkane fuel chemistry is 
needed. 

In the present work (ref. 4) experimental results on extinction of burning heptane pools were 
employed along with the asymptotic theory for the four-step description to infer what value of 7 
is needed to best fit the data. In this manner, it was found that 7 = 13 is close to the optimum 
selection. This suggests that revisions in equations (15) or (16) or in the propene and ethene 
consumption mechanisms are needed if agreement with extinction experiments is to be achieved. 

Computations were made not only with the four-step mechanism but also with a three-step 
mechanism, obtained from the four-step description by putting the H atom in steady state, and with 
a two-step mechanism, obtained from the three-step description by imposing partial equilibrium 
on step (8). There were large differences in the resulting predictions of heptane flame structures 
and extinction, much larger than for methane. For example, differences in the scalar dissipation at 
extinction exceeded a factor of ten, in contrast to being less than a factor of two for methane. It 
is concluded that the higher hydrocarbons are more sensitive to the kinetic descriptions and that 
therefore experiments on flame structures and extinction for these fuels may provide information 
of greater relevance to the chemical kinetics. 

Flame Structure and Extinction 

In the two-step approximation, the structure is as illustrated in figure 2, with the fuel-consumption 
zone thin compared with the oxygen-consumption zone, which in turn occupies an asymptotically 
small range of Z . In the three-step approximation, nonequilibrium of step (8) occurs in a portion of 
the oxygen- consumption zone adjacent to the fuel-consumption zone, provided that the Damkohler 
number for this step is large enough, which it is for the propene path but not for the ethene path. 
If the ethene path without modification is correct (which seems quite unlikely), then an alternative 
analysis, paralleling that of Bui-Pham, Seshadri, and Williams for the premixed methane flame 
with slow CO oxidation, will be needed, but it does not appear that the correct fuel chemistry for 
alkane diffusion flames will require this. In the four-step approximation, the thin fuel- consumption 
zone itself has substructure. There are two limiting behaviors of this substructure, as illustrated in 
figures 4 and 5. When H has a strong tendency to achieve steady state, as was found for methane, 
there is a thin radical-consumption layer within the fuel-consumption layer, as illustrated in figure 
4, where the nondimension al thickness of the fuel-consumption layer is and the parameter /i, (ref. 
4) involving ratios of reaction rates in the rate-ratio asymptotic analysis (RRA) is large. When this 
tendency is weak (/i small), the situation illustrated in figure 5 occurs, i n wh ich there is a diffusion 
flame between radicals and fuel within this inner layer — a diffusion flame within the diffusion flame. 
The results (ref. 4) show heptane to lie strongly in this last regime (^ -+ 0). 

Predictions of the location Zq and outer-solution temperature T° of the inner fuel-consumption 
zone are of interest. For heptane, contrary to the relative location illustrated in figure 2, the 
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two-step and three-step descriptions give Zq on the oxidizer side of Z c for most values of Xc> in 
particular at extinction. This disagrees with experiment and also with the predictions of the four- 
step mechanism, which place Z 0 on the fuel side, as illustrated in figure 2, even at extinction. The 
temperature T° of this fuel-consumption zone also is more realistic with the four-step mechanism, 
as illustrated in figure 6. Thus, in every respect, the four-step description is much better than 
simpler descriptions for higher hydrocarbons. 

From Xco equation (5) is used to calculate extinction diameters. Results are shown in figures 7 
and 8. Figure 7, for the propene route, clearly shows the large differences obtained with different 
degrees of reduction of the chemistry. The most reasonable values probably are close to the two- 
step prediction, fortuitously through cancellations of errors, but until heptane droplet extinction 
is observed experimentally, this conclusion is not certain. If the four-step predictions in figure 7 
are correct, then it is unlikely that droplet extinction can be observed experimentally, because the 
droplets will become too small. If 7 = 13 in equation (7), then the four-step predictions fall much 
closer to those shown for the two-step predictions. Comparison of the three-step results shown 
in figure 8 for the propene and ethene routes shows the extremely large difference in extinction 
diameters for these two descriptions (7 = 10 versus 7 = 2); calculations with two-step and four- 
step chemistry for the ethene route gave essentially zero extinction diameters. It seems likely 
that many more than two radicals are destroyed for each heptane molecule consumed, but further 
chemical studies axe needed to determine how many. 

Future Plans 


It is planned to investigate the fuel chemistry further, in an effort to determine better the values of 
the parameters a, /?, and 7. It is also planned to initiate similar asymptotic analyses for alcohol fuel 
droplets. Finally, microgravity experiments on alkane droplet extinction in different atmospheres 
are planned, to test predictions. 

References 


1. Choi, M.Y., Dryer, F.L., Card, J.M., Williams, F.A., Haggard, J.B., Jr., and Borowski, B.A., 

“Microgravity Combustion of Isolated n-Decane and n-Heptane Droplets,” AIAA Paper No. 
92-242, January, 1992. 

2. Card, J.M. and Williams, F.A., “Asymptotic Analysis of the Structure and Extinction of Spher- 

ically Symmetrical n-Heptane Diffusion Flames,” Combustion Science and Technology 84 
(1992)91-120. 

3. Card, J.M. and Williams, F.A., “Asymptotic Analysis with Reduced Chemistry for the Burning 

of n-Heptane Droplets,” Combustion and Flame, to appear, 1992. 

4. Card, J.M., “Asymptotic Analysis for the Burning of n-Heptane Droplets using a Four-Step 

Reduced Mechanism,” Combustion and Flame, submitted, 1992. 

5. Williams, F.A., “Combustion Theory,” Addison- Wesley, Menlo Park, CA, 1985. 

6. Kim, J.S. and Williams, F.A., “Structures of Flow and Mixture-Fraction Fields for Counterflow 

Diffusion Flames with Small Stoichiometric Mixture Fractions,” Combustion Science and 
Technology, submitted, 1992. 


287 



! 

j 

i 

I 

j 


i 


WVCCTIVt- REACTIVE- C0NVECT1VE- 

WfTUSIVC DIFFUSIVE DIFFUSIVE 

ZONE ZOMC ZONE 



o Zc i 


Figure 1.- Schematic diagram of the flame structure in the mixtare- 
traction coordinate far ooe-atep empirical Arrhenius chemistry, 
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Figure 2.- Schematic diagram of the flame structure 
in the mixture- fraction coordinate for two-step 
reduced chemistry. 
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Figure 3.- Schematic illustration 
of implications of theory concern- 
ing heptane flame structure in 
physical coordinates. 



Figure 4 .- Schematic illustration of inner-layer 
structure when there is a strong tendency for H to 
be in 6teady state. 



Figure 5.- Schematic illustration of inner-layer 
structure when there is a weak tendency for H to be 
in steady state. 
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temperature at extinction on the stoichiometric 
mixture fraction (varied by dilution) for the propene 
path at normal atmospheric pressure and 
temperature. 


EXTINCTION DIAMETER ESTIMATES 



extinction diameter on the oxygen mole fraction in 
the ambient gas for the propene path. 


EXTINCTION DIAMETER ESTIMATES 



Figure 8.- Predicted dependence of the heptane 
extinction diameter on the oxygen mole fraction in 
the ambient gas for both propene and ethene routes. 
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COMMENTS 


Question (Ray Edelman, Rockwell International): I would expect burning rale and extinction, in particular, 
would depend strongly on soot formation and the attendant radiation loss on reducing temper, ilure. This 
might be enhanced at high pressures. Could you comment on what your data shows and how you reconcile 
this with your analyses? Do competing effects balance this with enhanced kinetics at increasing pressure? 

Answer : Our studies are at normal atmospheric pressures and at reduced pressures, not at high pressures. 

For our small single droplets under these conditions, radiant transfer associated with soot formation probably 
is negligible. At high pressures it will become more important, but temperature reduction through radiant 
energy loss still could be a perturbation because of the small size. 
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